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In the SUSY SO(10) GUT context, we study the exclusive processes B — > K^l + l~(l = fJ,,r). 
Using the Wilson coefficients of relevant operators including the new operators Q-[\ which are in- 
duced by neutral Higgs boson (NHB) penguins, we evaluate some possible observables associated 
with these processes like, the invariant mass spectrum (IMS), lepton pair forward backward asym- 
metry (FBA), lepton polarization asymmetries etc. In the model the contributions from Wilson 
coefficients C'q 1 2 , among new contributions, are dominant. Our results show that the NHB effects 

are sensitive to the FBA, dL/ds, and dT/ds of B — * K^*'t + t~ decay, which are expected to be 
measured in B factories, the deviation of dT/ds in B — » Kfi + fi~ can reach 0.1 from SM, which could 
be seen in B factories, and the average of the normal polarization dN/ds can reach several percent 
for B — > Kfi + /j,~ and it is 0.05 or so for B — > Kt + t~ , which could be measured in the future super 
B factories and provide a useful information to probe new physics and discriminate different models. 



I. INTRODUCTION 



The rich flavor changing neutral current processes B — > K^*H + l~ have been the sharper focus since these decays 
are potential testing grounds for the SM at loop level and hoped to probe the new physics beyond the SM. Recently 
exclusive measurements have been done by Belle and BaBar and the following results for the branching ratios of the 
B -> Ki + £~ and B -> K*£+£- (I = e,n) decays are announced 0,0: 



Br(B -> Kl + t 



(4.8±J;g ± 0.3 ± 0.1) x 10~ 7 [Belh 



(0.651^3 ± 0.04) x 10~ 6 [BaBar] 



(1) 



Br{B -> K*£ + £-) = 



[ll. 5±|;| ± 0.8 ± 0.2) x 10~ 7 [Belle] 



(0.88±° 



0.33^i 
29 J 



X 10 



-6 



[BaBc 



which imply 



Br(B -> Kl+l-) wmld ^ vc = (0.54 ±0.09) x lO" 6 
Br(B K*l+l-) woM ^ vc = (1.04 ±0.22) x lO" 6 



(2) 



In addition, the FBA for B — > K*l + l~ have firstly observed at Belle collider 0. 

The B — > K(*h + l~ decays, induced by b — > sl + l~ transition at quark level, are experimentally easier to measure 
than inclusive processes B — > X s l + l~ . From the theoretical point of view there are large uncertainties, which come 
mostly from the decay form factors, to make predictions for the exclusive processes. At present, the knowledge of 
form factors lacks a precisely non-perturbative solution. A number of papers are dedicated to calculating the form 
factors with various appropriate methods [USUI El Among them, the QCD sum rules on the light-cone (LCSRs), 
which deals with form factors at small values of s, the momentum transfer to leptons, is a complementary to lattice 
approach and has the consistence with perturbative QCD and the heavy quark limit. In this paper, we will use the 
form factors calculated by the LCSRs [Tj|. 

The measurement of invariant mass spectrum, forward-backward asymmetry, and lepton polarizations are efficient 
tools to establish the new physics. There are a great deal of studies for processes B — > K^l + l~ in theory. A 
model independent analysis has been carried out in J?] , Q and a lot of pap ers per form the investigation in many new 
physics scenarios [SElllEEllil J HJliJ O llSrand some works [11 IHlHj are dedicated to the double lepton 
polarization. It has been pointed out [Tj, [HJ Eg that in the some types of two-Higgs-double model and SUSY models 
the neutral Higgs bosons have sizable contributions to these decays (for I = fj,,r) at large tan/3. In Ali et.al 
calculate these quantities in five scenarios of supersymmetric model with assuming no additional phase. Kruger's 
studies 0| are focus on the CP violation of B — * Kl + l~ in the model with additional CP phases and an extended 
operator basis. In ref. (I3j only the Higgs pen guin s with chargino-stop propagated in the loop are considered. 

The rapid progress in neutrino experiment [22j requires new physics to provide a theoretical explanation. Moti- 
vated by neutrino observations, a number of SUSY SO(10) models have been proposed p3l I24I EH |2(| and some 
phenomenological consequences of the models have been discussed [23|, |23, 123, |2!| y(j In SUSY SO(10) GUT models, 
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there is a complex flavor non-diagonal down-type squark mass matrix element of 2nd and 3rd generations of order 
one in the RR sector (i.e., 8$3 R non-zero with tffiP* = (^Jfljj) 23 / 771 !' wnere (^^ RR )ij 1S the flavor non-diagonal 
squared right-handed down squark mass matrix element and rriq is the average right-handed down-type squark mass) 
at the GUT scale [25| which can induce large flavor off-diagonal couplings such as the coupling of gluino to the quark 
and squark which belong to different generations. These couplings are in general complex and may contribute to the 
process of flavor changing neutral currents(FCNC). For specific, we use the SUSY SO(fO) model described in ref. [25| . 
The details and a simple description of this model can be found in Refs. [2f|,|3(|. In this paper, we investigate exclusive 
decay B —> K^l + l (I = M,r) in the context of SUSY SO(IO) GUT. It is well-known that the effects of the counter- 
parts of usual chromo-magnetic and electro-magnetic dipole moment operators as well as semileptonic operators with 
opposite chirality are suppressed by m s /mb and consequently negligible in SM. However, in SUSY SO(10) GUTs their 
effects can be significant, since S^ 11 can be as large as 0.5 [25J. Furthermore, 5^ R can induce new operators, the 
counterparts of usual scalar operators (Qi,2, for their definitions, see below) in SUSY models, due to NHB penguins 
with gluino-down type squark propagated in the loop. We include the contributions of these counterpart operators 
and find that indeed they are dominant in the SUSY SO(I0), using the MIA with double insertions to calculate 
Wilson coefficients of operators. The aim of our paper is make an analysis of the SUSY contributions, in particular, 
the contributions of neutral Higgs bosons, to the exclusive decay B — > K^l + l^(l — /i, r) in the context of SUSY 
SO(10) GUT. 

The paper is organized as follows. In section 2, we present the effective Hamiltonian and hadronic matrix elements 
of relevant operators in terms of form factors. In section 3, the expressions of observables are given. In section 4, we 
give the sparticle mass spectrum using the revised ISAJET. We make the numerical analysis and draw the conclusion 
in section 5. 



II. EFFECTIVE HAMILTONIAN AND FORM FACTORS 



In the SUSY SO(I0) GUT, after integrating the heavy degree of freedom from the full theory, the general effective 
Hamiltonian for b — > sl + l~ can be written as follows: 



'75 



v tb v t * 



i=3 



i=l 



(3) 



where Oi(n)(i = 1, • ■ • , 10) are dimension-six operators and Cj(/z) are the corresponding Wilson coefficients at the 
scale \x [3l|. The additional operators Qi{i = I, - • - ,8) come from the neutral Higgs exchange diagrams and their 



definitions are given as [la, |3 



Qi 
Q-2 

h(4) 



16n 2 

e 2 
16n 2 

-,2 



(s a L b a R )(ll) 
(s a L b R )(h 5 l) 



L)> 



05(6) = if^(*z&£)(£ 



'5(6) 
Ql 



16n 



16n- 



■(■^b R )(J2f L ^,q a R ) 



(4) 



and the corresponding Wilson coefficients can be found in |33j . The primed operators, the counterpart of the 
unprimed operators, are obtained by replacing the chiralities in the corresponding unprimed operators with opposite 
ones. The explicit expressions of the operators governing B — > K^*H + l" are given as: 



7 
Og 



16tt 2 

e 2 

f 67T 



m b {sa^P R b)F^, 
■{s ltl P L b){hH), 



0' 7 = j^m b (sa^P L b)F^, 



OL 



16tt 2 



(s^P R b)(h»l), 



3 



2 2 

Oio = i^(57/ift6)(iy7s0, O' 10 = I £- 5 (a 7M P fl 6)(zV 75 0, 

2 2 

2 2 

Q 2 = ^(sPfitJ^l), Q2 = ^(sft&Xi-ysO (5) 

where Pl.jj = (1 T 7s)/2. From the above Hamiltonian, we get the decay amplitude of 6 — > sl + l~: 

Gfo> 



M(b sl+r) = - -^VtbV t * s [cf f [s^Lb][l^l] + Cio[s7^HV75*] 



-2m b C 7 e// [sia^Ci?6][I 7 ^7] + C 0l [5i»](H] + Cq 2 [I 75 /] + (Ci(m 6 ) «-> C{(mt)) \ (6) 



where s = q 2 ,s = -V, q = p B ~ Pk(*"> ^ s the momentum transfer. The Wilson coefficient Cg (/z) and C 7 are 
defined as: 

C 9 e// (A* ; s) = C 9 (/i) + F(/x, 3) + ^C(M)S Vi=v , (ls) ,,, v , (6s) fc i ^^^| j i!^Tr v . ' ^ 

C 7 e// = C r -C 6 /3-C t (8) 

where Cg^(/z) contains the long distance effects associated with real cc in the intermediate states B — » KJ/ip(ip') — > 
Kl + l~, which can be expressed as the last term in Eq.J7J), as well as the short distance contributions. The function 
Y([a, s) comes from the one-loop contributions of the four-quark operators and its explicit expression can be found 

in |34j |. The C 9 e ^(/z) and C 7 e ^ can be obtained by replacing the unprimed Wilson coefficients with the corresponding 
primed ones in the above formula. 

In virtue of the form factors in [l4|, the hadronic matrix elements in the B — > Kl + l~ decay can be expressed as: 

< K(p)\s ltl b\B(p b ) > = f+(a) Pll + (9) 

< K(p)\S* MV ?(l + To)b\B( Pb ) > = 1 {p^s - q^ml - m 2 K )} h [ s) (10) 

m B + m K 

Using equations of motion, we obtain 

< K( PK )\sb\B( PB ) > = f (s) . i I I ! 

m s - m b 

For B — > K*l + l~ , the form factors are defined as follows. 



<K*(p K ,,e)\s^(l±'j 5 )b\B(p B ) > = e e *>>^,— — ±ie*Jm B +m K *)A 1 (s) 

m B + m K . F 



. / * v Ms) 

T»P M ( £ Pb)- 



m B + m K * 



Ti%{e*p B ) 2 ^{A z {s) - A (s)), (12) 



< K*(p K .,E)\sa^q v (l± l5 )b\B(p B ) > = ie ll „ IMT e" , p p B p' I K .2T 1 (s) ± e;T 2 (s){m B - m 2 K 



T(e>fe T 2 (s)+T 3 (s)— — ± (e*p B ) qfl T 3 (s) (13) 



and 



m B - m K > 



<K*(p K .,e)\s(l± l5 )b\B(p B )> = Ti(e* PB ) 2m *' A (s) (14) 

mb + m s 



by means of equations of motion. 

The q 2 dependence of the form factors can be parameterized as 

F(s) = F{0)exp{ds + c 2 s 2 + c 3 s 3 ) 

where related parameters are given in the Table. 4 of [l4j . 
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III. THE FORMULA FOR OBSERVABLES 



From Ea. l|3ll4|l . we can write the decay matrix elements as 

Gfoi 



A 



where for B — ► Kl + l decay: 



V tb V t *m B [T^Q-fl) + TlQrisl) + SQl)} 



S = S x {s) 



and for B — > K*l + l decay: 



T l = Ms) e wctt3£* p PBPK* ~ iB ( g yt + iC{s)(e* ■ p B )p^, 

Tl = E(s)e„ pa pe*?p B f K ,-iF(s)e*+iG(s)(e* ■ p B )p^ + iH{s){e* ■ p B )q^ 



S = i2m K *(e* ■ p B )S 2 (s) 
where p = ps + Pk(*) > 5 = Vb — Pk(*) j ^ = ""^S P — an d the auxiliary functions are defined as: 



A'(s 
C'(s 
D'(s 

Si (« 
A(s 
B(s 
C(s 

E(s 
F(s 
G(s 

H(s 

S 2 (s 



1 + ttik 

{C 10 + C' w )f + (s), 
(C 10 + C( )/_(s) 



1 — rh 2 K 



2rhi{rhi, — rh s ) 



(C Q2 +C' Q2 )f (s), 



1 + niK* s 

(1 + m K ,)[Cl ff (s) - C' 9 eff (S)]A X (S) + ^(1 - m 2 K ,)(C e 7 ff - C 7 eff )T 2 (s), 



(l + m*.) 1 "" 
2V(s 



1 + rh K 
2rriK* 



(Cio + C 10 ), 

)(Cio - Cio) 
(C 10 - C' w )A 2 (s), 



1 + 

(1 +m K .)(C 1Q -C' w )A 1 (J), 
1 



(C 10 -<7i )(A 3 (s)-A (s)) + 

s 

^^^)(C* Ql -C Ql ), 



m ; (m fc + m s ) 



(Cq 2 -C' Q2 )A (s) 



(15) 



(16) 



(17) 

(18) 
(19) 
(20) 

(21) 

(22) 

(23) 

(24) 

(25) 
(26) 
(27) 

(28) 

(29) 



where /_(s) = mB ~ mK (fo(s) - /+(«)), A 3 (s) = m f+". K * Ai(a) - m f A 2 (s). The above results reduce to those in 
ref. ^3 ^ au Cj' = 0, as expected. It is worth to note that the final term in cq. (|15f) vanishes if one does not include 
the NHB contributions. 



A. The dilepton invarient mass spectra and differential FBA 

The kinematic variables s, u are defined as: 



s = q 2 = {p++p-f, 

U = (PB - P-f - (pB - P+f 



(30) 
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Here we choose the center of mass frame of the dileptons as the frame of reference, in which the leptons move back 
to back, and the momentum of B meson makes an angle 9 with that of l + . u can be written in terms of 9: 

u = —u(s)-cos9 = —ii(s)z,z = cos9, 
u(s) = 



m 



KM 



A = 1 

The phase space is defined in terms of s and z: 

(2rhi) 2 < s < (1 - m KM ) 2 



A(l -4— l -),V = \/l -4- 
s V s 

s 2 - 2s- 2m 2 w (l + s) 



?,4 



-1< z < 1 



(31) 



(32) 



Keeping the lepton mass and integrating over u in the kinematic region, we can get the dilepton invariant mass 
spectra (IMS): 



dT K ( K ' 
ds 



G > 2m \ Vth V t *fu{s)D K ^'), 



D K = (\A 



2 10 7T 5 
l\2 i 



Cf) { X-^l ) + ]Sl \ 2 (s-4mf) 



+ \C , \ 2 4mf(2 + 2rh 2 K - s) + Re(C D'"<)8m 2 (l - rh 2 K ) + \D'\ 2 Am 



Is. 



(33) 



(34) 



D K * = 



I4! SA (1 + 2 % + &su(s) 2 + \S 2 \ 2 (s - 4m 2 )A 
3 s 3 



Arh 2 K , 
X 

4rh 2 K , 



\B\ 2 (X- 



|C| 2 (A- 



ii{s) 2 



+ 8m 2 K ,{s + 2m 2 )) + \F\ 2 (X- 



u{s) 2 



+ 8m 2 K ,(s - 4to 2 )) 



) + \G\ 2 X- 



1 



2rh 2 K , 



Re(BC^)(l 



+Re(FG^) ( (1 - 



s){\- 



s)(X- 
u(s) 2 



u(s) 2 
3 

u(s) 2 



+ 4m 2 (2 + 2m 2 K , - s) 



4to 2 A 



rr> 



K* 



The differential FBA is defined as: 



-A [Re(FH^) - Re(GH^){\ - rh 2 K ,)} + \H\ 



K* 



A 



FB(s) 



JO dsdu J —u{s) dsdu 

JO dsdu J —u(s) dsdu 



(35) 



(36) 



According to the definition, it is straightforward to obtain the expressions of FBA in the exclusive decays: 
• B -► Kl+l- 

^fB D K = -2m l &(8)Re(S l A'*) 
ds 



(37) 



• B -> K*l+l~ 



rlA K * 
a/i FB jjK" 

ds 



(s) <^ s[Re{BE^) + Re(AF^)} 



mi 

rh K * 



[Re(S 2 Bl)(l — s — m 2 K .) - Re(S 2 &)\] 



(38) 



Seen from the Eqs. {SU and (JSHJ, the functions D'(s), 5i(s), H(s), and S 2 (s), which come from the 

contribution of NHBs, enter the IMS and FBA. Hence, the effects of NHBs will manifest themselves in the numerical 
results of these formula. In particular, Eq. I|37[) shows that the FBA in B — > Kl + l~ vanishes if there is no NHB 
contributions and from Eq. (|38|l it follows that the NHB contributions change the position of the zero-point of the 
FBA in B — > K*l + l~ . As pointed out in ref. [20|. in an untagged sample, the FB asymmetry for unpolarized leptons 
vanishes. Once the flavor of the decaying b-quark is tagged, one can measure the unpolarized FB asymmetry which 
is an important observable to discriminate new physics from the SM, as we noted above. 
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B. The lepton polarization 

In this subsection, we will present the analytical expressions of lepton polarization. We define the three orthogonal 
unit vectors in the center mass frame of dilepton as 

<?L = P+, 

Pk x p+ 



e N 



\PK XP+\' 

<§t = e N x e L , (39) 

which are related to the spin of lepton by a Lorentz boost. Then, the decay width of the B — > K (*H + l~ decay for 
any spin direction h of the lepton, where h is a unit vector in the dilepton center mass frame, can be written as: 

= ^ (§) I 1 + (P^l + P N e N + P T e T ) • A] (40) 

where the subscript "0" denotes the unpolarized decay width, Pj, and Pr are the longitudinal and transverse polar- 
ization asymmetries in the decay plane respectively, and Pjv is the normal polarization asymmetry in the direction 
perpendicular to the decay plane. 

The lepton polarization asymmetry Pi can be obtained by calculating 

_ dT(h = e t )/ds - dT(n = -e t )/ds 
dT(h = ei)/ds + dT{n = —ei)/ds 

By a straightforward calculation, we get 
• B -> Kl+l- 

P*D K = ^V{\Re(A'C n )-3mi(l-m 2 K )Re{C'' i S 1 )-3m l sRe(D'' i S 1 )} , (42) 



P%D K = 
n 2 



{-Im(A'Sl) + 2m;Im(C"£>' t )} , (43) 



P«D K = jm, [(1 - m 2 K )Re{A'C*) + §Re(A'D*)] + ( * j™H R e {C'S\)} (44) 

B -> K*l+l 



PfD K ' = v\—Re(AE^+ {X + 12 f K ' S) Re(BF^ 



Mm|^i) fe(BGt+CFt)+ _£_ fe(CGt) 

2m;A 



[Pe(P £|) - SRe(HSl) - (1 - m^)^(G5 2 f )] J , (45) 

Attik* [_ttik* 
+Pm(PP t )(l - m 2 K , - s) - Im{GH^)X] 
+2m K ,rhi[Im(BE j! ) + Im(AF r )} 

-(1 - m 2 K . - s)Im(BSl) + \Im(CSl)\ , (46) 

P f D K- = Z^{4JPe(^t ) 
4Vs 

+ (1 - m f - ~ s) [-Pe(PPt) + (1 - mi.)-Re(PGt) + jp e (ppt)l 

+ ^-[Pe(CP t ) - (1 - mi,)Pe(CG t ) - sPe(CP"+)l 
m 2 K , 

+ (s_4m i ] _ ^ _ S , R , FS U _ X Re(GSt)]\ (47) 
m K *mi J 

One can see from Eq. (|43|l that P/v = for the decay B — > Kl + l~ in the SM because C{ = in the approximation 
of m s /mi, — 0, Cq| 2 = 0, and C\q is real in the SM. Thus, a non zero normal polarization asymmetry in B — > Kl + l~ 
would signal the existence of new physics. 
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FIG. 1: The mass spectrum versus mo for fixed Gev, fem/3=40, 523fli{ = (0.04-0.03i), and sign(pi)=+l without the 

constraints from the low energy experiments imposed, (a) is for Ao = 0. (b) is for ^4o=-1000. 



IV. MASS SPECTRA AND THE PERMITTED PARAMETER SPACE 



To see the impact of the induced off-diagonal elements in the mass matrix of the right-handed down-type squarks 
on B rare decays and simplify the analysis, we assume that at the GUT scale (Mg) all sfermion mass matrices except 
the right-handed down-type squark mass matrix are flavor diagonal and all diagonal elements are approximately 
universal and equal to jtiq. The 2-3 matrix element of right-handed down- type squark mass matrix is parameterized 

by 5H R EE Which can be treated as a free parameter of order one. Furthermore, we have a universal gaugino 

mass Mi/2, a universal trilinear coupling Ao and a universal bilinear coupling Bo at Mq. Taking account to the 
radiative electro-weak (EW) symmetry breaking, finally we have five parameters (mo, M\/2, Ao, S^* 1 , tan/3) plus the 
sign of fx as the initial conditions for solving the renormalization group equations (RGEs). 

We require the lightest neutralino to be the lightest supersymmetric particle (LSP) and use several experimental 
limits to constraint the parameter space, including l)the width of the decay Z — > XiXi 1S l ess than 4.3 MeV, and 
branching ratios of Z — > X1X2 anc ' Z — > X2X2 are l ess than 1 x 10~ 5 , where \i is the lightest neutralino and X2 * s 
the other neutralino; 2) the mass of light neutral Higgs can not be lower than 111 GeV as the present experments 
required; 3) the mass of lighter chargino must be larger than 94 GeV as given by the Particle Data Group jljj; 4) 
sneutrinos are larger than 94GeV; 5) seletrons are larger than 73GeV; 6) smuons larger than 94 GeV; 7) staus larger 
than 81.9 GeV. 

In the numerical calculation, we use the revised IS A JET. We find that the parameter #33 d° es not receive any 
significant correction and the diagonal entries of mass matrices are significantly corrected, which is in agreement with 
the results in Ref. |4l|. We scan m n , M X j 2 in the range (100, 800) GeV for given values of A , tan (3 and sign(/i)=+l*, 
with the constraints from the relevant low energy experiments such as B — > X s j, B s — > fi + fi~ , etc. (for the detailed 
discussions of constraints, see VB). 

For an illustration, we present the mass spectra without and with the constraints from the low energy experiments 
in Figs. 1 and 2, respectively, where (a) and (b) are for Aq — 0, — lOOOGeV respectively. One can see from the Figs. 
^ 121 that the mass spectrum lifts when Ao increases and when the constraints from the low energy experiments are 
imposed the masses of sparticles are larger than those without the constraints, as expected. 



V. NUMERICAL ANALYSIS 



In this section, we will discuss the numerical results and make an analysis. 



* In the case of sign(/^)= -1, the constraint from B —> X s "/ on the parameter space is too stringent, in particular, for large tan /3 1321144.1451 
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FIG. 2: The mass spectrum versus mo for fixed 71^1/2=600 Gev, tan/3=40, <523ii.R = (0.04-0.03i), and sign(/x)=+l with the 
constraints from the low energy experiments imposed, (a) is for Aq — 0. (b) is for J 4o=-1000. 



Parameters in our calculation are list 



A. Parameters input 



m b = 4.8GeV, m c = lAGeV, 
M B = 5.28GeV,M J/4 , = 3.10Gey, 
T B = 4.22 x 10 

r(j/v>->i + r) 



13 GeV, 

= 5.26 x 10- 6 GeV, 



m s = 0.2GeV,m p = 0.1057GeT/, 
M$> = 3.69GeV, M K , = 0.89GeV, 
Tj/^ = 8.70 x 10~ 5 GeV, 
T(ip' -> l+r) = 2.14 x 10~ 6 GeV 



m T = 1.7769GeV, 

M K = 0.49Gey, 

I> = 27.70 x 10~ 5 GeF 



B. Constraints from experiments 

In our calculation, we consider the constraints from B — > X s j, B s — > AM S , £> — > K(*H + 1~, and r — > /x 7 . 

The leading order _B S — > X s ~/ branching ratio normalized to Br(B — > X c eD) can be written as 

Br(B -> X s7 ) = ^pL|^|3 B r(B - X c eV){\C 7 {m b )\ 2 + |^(m b )| 2 ) 
717(2) V cb 

where VI = mP° le /m p b ° le , f(z) is the phase space function. We take 2 x 10~ 4 < Br(B -> X s7 ) < 4.5 x 10~ 4 , 
considering the theoretical uncertainties. The B — ► X s7 make a direct constraint on \Cr(m b )\ 2 + \C 7 (m b )\ 2 . The 
single insertion term in C 7 (mb) are more severely constrained than 8%£^ R , due to the strong enhance- 

ment factor rrig /mi, associated with single J^^^^ insertion term in C 7 \mh). Because the double insertion term 
^dLL(RR) gdLR(LR*) ^ enhanced by m§ /m&, #23^ is constrained to be order of 10~ 2 if the left-right mixing 
of scalar bottom quark 5^q R ^ RL ' is large (~ 0.5). Nevertheless, in the large tan/3 case the chargino contribution can 
destructively interfere with the SM (plus the charged Higgs) contribution so that the constraint can be easily satisfied. 
The branching ratio Br(B — > is given as 



Bi(B a ^»+iT) = ^f-™| 3 r B J 2 JA^yi^^[(l-4m 2 )|G Ql (m,)-G^(m b )| 2 + 

|Gq 2 (to 6 ) - C' Q2 (m b ) + 2m(C w (m b ) - C' w (m b ))\ 2 } (48) 

where fh = m^/rriB,. With large Gg| 2 , Br(B — ► /i + /i~) can have large enhancements [3(|. The new Dq experimental 



upper bound of Br(B s — > /z+/z ) is 4.6 x 10 7 [33 at 90% confidence level. It gives a stringent constraint on GqJ 
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TABLE I: The Wilson coefficients for the two cases in the SUSY S'O(IO). The SM values also are listed for comparison. The 
values in bracket are for I = r. 

Aq Cqi Cqi ( 2Sl ^hi 

SM 

Ap = 0.074+0.1(1.252+0.0011) -0.013+0.008I(-0.213+0.128I) -0.075+0.I(-1.267-0.001I) -0.013+0.008I(-0.216+0.129I) 
A = -1000 0.106+0.1(1.775+0.0021) -0. 247+0. 242I(-4. 148+4.0741) -0.107+0.I(-1. 797-0. 0021) -0.250+0.246I(-4.202+4.128I) 

Aq Cy C 7 B ^ Co, Cg Cio C'iq 

SM -0.313 +4.344 -4.669 

A = -0.225-0.1 -0.02-0.011 4.277+0.1 -0.-0.002I -4.717-0.1 -0.001+0.0191 
A = -1000 -0.219+0.1 0.039-0.0381 4.275+0.1 0.011+0.0721 -4.732-0.1 -0.075-0.6701 



and consequently the parameter space of the model. At the same time we require that predicted branching ratios of 
B — > X s fi + /i~ and B — > falls within 1 a experimental bounds. 

We also impose the current experimental lower bound AM S > 14.4PS" 1 The 5^ R contribution to AM S is 
small because it is constrained to be order of 10~ 2 by Br(B — > X s j). The dominant contribution to AM S comes from 
5^ R insertion with both constructive and destructive effects compared with the SM contribution, where the too 
large destructive effect is ruled out, because SM prediction is only slightly above the present experiment lower bound. 

Furthermore, as analyzed in ref. |4l| . there is the correlation between flavor changing squark and slepton mass 
insertions in SUSY GUTs. The correlation leads to a bound on S^ 1 * from the rare decay r — > fij. We update the 
analysis with latest BELLE upper bound of Br(r -> /ry)< 3.1 x 1(T 7 at 90% confidence level in the SUSY SO(10) 
model. 



C. The Numerical results and Conclusions 

We will focus on the parameter space at large tan(3. The reason for this is that in the large tan/3 region of 
parameter space the contributions of NHB exchange become very important for quark level semi-leptonic transitions 
b — ► sl + l~ when the final state lepton is either a muon or tau |32l l45j . In numerical calculations, we take tan/3=40, 
sign(/i)=+l, and Aq = 0, —1000 and get the sparticle mass spectrum and mixing at the EW scale. Using the resulted 
Wilson coefficients, we calculate the IMS, FBA and polarization asymmetries of processes B -► K^l+l- under the 
constraints on S^ 11 from the all relevant experiments as discussed in subsection VB, whose phase varies from to 2tt. 

The Wilson coefficients Cq| and Cq^ come from NHB exchanging. Specially, we are interested in the case of 

maximal enhancements of Cq| 2 „ Through scanning the parameter space under constraints, we find, for Aq = 0, when 

m =800GeV, M 1/2 =400GeV (for A =-1000, M 1/2 =500GeV), C%1 , have maximal values. The obtained C%\ , and 
other relevant Wilson coefficients in the two cases and in SM are listed in Table I. From Table I, we can know: (i) 
the NHB contributions in the case of Aq = —1000 are larger than those in the case of Aq — except for Cg; (ii) 
the Wilson coefficients C'q 1 2 of primed operators are dominant, which is due to the presence of #33 of order one 
at the high scale in the SUSY SO(10) model, and their imaginary parts are sizable, which contribute to the normal 
polarization, a T violating observable. 

The figures for the dependence of observable on s with and without long-distance contributions are presented in 
Fig. (|3l7|) in the case of Aq = —1000, where the solid lines denote the all contributions (W ± ,H ± , chargino, gluino, 
neutrilino propagated in the loop) including the NHB contributions; the dot lines present the SM contribution plus 
only the NHB contributions, and the dot-dashed lines are for the SM contribution. We also calculated the dependence 
of observable on s in the case of Aq = and give the results for Aq = when the two cases have a sizable difference. 

The IMS of the process B — > K^l + l~ is given in FigEI where the left two figures are for B — > K^fi + n~ and 
the right for B — > K^*\ + t~ . For the case of B — > Kfi + fi~, we can see that, at the low s region, there are some 
enhancement from NHB contributions. Compared to the decay B — > if/i+u - , the IMS of B — > K*fi + fi~ deviates 
from the SM prediction sizably in the whole region of s. Nevertheless, the SUSY effects are small compared with the 
SM for the decay B -► K^t+t~. 

The Fig. His for the FBA (dA/ds) of the decay B -> K^l+T, where the left two fi gures are for I = \x and the 
right for I = r, like that in FigEl As it is known, the FBA (dA/ds) of B -> Kl+l~ in the SM is zero. Since FBA 
arises in the SUSY models only when NHB effects are taken into account, it provides a good probe to test these 
effects. Our numerical results show that the average of FBA in B — > K/i + fi~ can reach only 0.001 which is too small 
to be observed. The average of FBA in B — > Kt + t~ can reach -0.1 and 0.05 for the case of Aq = and Aq = —1000, 
respectively. (The reason why FBA in the two cases has the opposite sign is that the sign of function Si in these 
two cases is opposite.) So, the 10 10 — 10 11 Bd Bd pairs per year, which is in the designed range in the future super 
B factors with B hadrons 10 10 — 10 12 per year [43, is needed in order to observe the FBA with good accuracy. Our 
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results show that the SUSY effects show up at the low s region for the FBA of B — > if*/i + /i~ and the deviation from 
SM is 0.05 or so. It is worth to note that there is a sizable change of the position of the zero-point of the FBA in 
B — ► A*/i + /i~ in the SUSY SO(10) model, as it can be seen in Fig. 5, which could be tested in the future experiments 
with high precision. For FBA in B — > K*t + t~, the deviation from the SM is about several percent. The average 
of FBA of B — ► K*t + t~ can reach 0.3 in the case of A a = 0. To observe the FBA in B — ► K*t + t~ decay at la 
level, the required number of events is 1.1 x 10 8 . The number of BB pairs that is expected to be produced at B 
factories is about N ~ 5 x 10 8 . Therefore the FBA in B — > K*t + t~ could be observed at B factories. Hence, with 
the enhancement of experimental precision and statistics, the measurements of FBA would provide more data and 
effectively pin to the NP effects. 

Now, we turn to discuss the lepton polarization. We present the longitudinal, transverse and normal polarization 
in the Fig. {51517) 1 for B -> K^l + l~ decay. As it can be seen in Fig. [5Jand|5J the dL/ds of B -> K^ti + fi~ is not 
sensitive to the NHB effects, while for dT/ds of B — > K{K*)fi + \i~ , the deviation from SM can reach 0.1 (0.05). As it is 
expected, the contribution from the t + t~ channel is much larger than that from the one. For B -> K^t+t-, 

the NHB contributions are manifest and dominant, and both dL/ds and dT/ds are significantly different from SM. 
And the dL/ds of B — > At + t~ can even reach 0.6. Thus, the NHB effects are sensitive to B — > if (*)t + t~ and will 
be observable at B factories. 

The dN/ds of S — > K^-*'l + l~ decay are given in Fig. UJ The average of dN/ds can reach several percent for 
B — > K/i + fi~ which could be observed in the future super B factories, while it is the order of 10~ 3 for B — ► A*/i + /i~ 
which can not be observed even in designed super B factories. The average of dN/ds in _B — * Kt + t~ is 0.05 or so. 
For i? — > K*t + t~ , the deviation from SM is a few percent. As noted above, the Wilson coefficient Cio is real and 
C' w , Cq. = (precisely speaking, they are negligibly small) in the SM so that dN/ds = in B — > K(J, + fJ,~~ in the SM. It 
is still true in the minimal super gravity model (mSUGRA) and SUSY models with real universal boundary conditions 
at the high scale 01 . In the SUSY SO(10) model we considered, the complex flavor non-diagonal down-type squark 
mass matrix element of 2nd and 3rd generations of order one at the GUT scale induces the complex couplings which 
lead to the complex Wilson coefficients and consequently the non zero normal polarization of B — ► K/i + fi~ . Therefore, 
the measurements of the CP violating (as usual, the CPT invariance is assumed in the paper) normal polarization 
in B — > Kl + l~ could discriminate the SUSY SO(10) model (and other SUSY models with the flavor non diagonal 
complex couplings) from the SM and mSUGRA. 

In summary, we have carried out a study of SUSY effects, in particular, the neutral Higgs bosons contributions to 
the IMS, FBA and polarization, in the exclusive decay B — > K^l + l~ (I = fi,r) in the SUSY SO(10) model, taking 
account of the constraints from existing experimental data such as b — ■» 57, AM S , Br(B — ► K (*)/i + /i~), r — > /17 as 
well as the upper bound of Br(B s — > /i + /i~). Our main findings can be summarized as follows: 

• The IMS of the process B — > K^*'fJ, + fi~ can sizably deviate from the SM. 

• The FBA comes only from NHB contributions in B — > Kl + l~ and its average for I = fi is nonzero but too small 
to be ovserved. However for B — ► Kt + t~ , it is the order of 10%, which should be within the luminosity reach 
of coming B factories. The SUSY effects show up at the low s region for the FBA of B — » K and the 
deviation from SM is 0.05 or so. Moreover, there is a sizable change of the position of the zero-point of the FBA 
in B — > K*fi + fi~ , which can be used to discriminate the model from the SM. 

• The average of dN/ds can reach several percent for B — * A"/i + /i _ and it is 0.05 or so for B — * Kt + t~, which 
could be measured in the future super B factories and provide a useful information to probe new physics and 
discriminate different models. 

• The longitudinal polarization, dL/ds, of B — > K^*' fi + fj," is not sensitive to the NHB effects. However, for the 
transverse polarization, dT/ds, of B — > A(AT *)/i + /i - , the deviation from SM can reach 0.1 (0.05) which could 
be seen in B factories. For £? — ► A(*'t + t~, the NHB contributions are manifest and dominant, and both dL/ds 
and dT/ds are significantly different from SM. And the dL/ds of B — > Kt + t~ can even reach 0.6, which can 
be measured in B factories. 

Therefore, the experimental investigation of observable, in particular, FBA and the polarization components, in the 
B — > K (*>l + l~ decays in the present B factories and future super B factories can be used to search for SUSY effects, 
in particular, NHB effects, in SUSY grand unification models. 
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FIG. 3: The IMS of the process B -> K ( *H + r for A = -1000. The solid line(black), dot line(green), and dashed-dot line(blue) 
represent the all contributions included, the SM contributions plus only the NHB contributions, and the SM contributions, 
respectively. Both the total (SD+LD) and the pure SD contributions are shown in order to compare. In the fig. we write "s" 
in stead of " s" for simplicity. 
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FIG. 4: The FBA of the process B — ► K^*H + l for Ao = —1000. The line conventions are the same as those in Fig. 
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FIG. 5: The dL/ds of the process B -> K M l + r for A = -1000. The line conventions are the same as those in Fig. |31 




FIG. 6: The dT/ds of the process B -> K { *H + r for A = -1000. The line conventions are the same as those in Fig. ED 




FIG. 7: The dN/ds of the process B -> K^l+T for A = 



— 1000. The line conventions are the same as those in Fig. [3] 



